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The Interference between Virtual Photon and 1 Charmonium in e + e~ 
Experiment 

P. Wang a *, C. Z. Yuan a , X. H. Mo a ' b , D. H. Zhang a 

institute of High Energy Physics, P.O.Box 918, Beijing 100039, China 

b China Center of Advanced Science and Technology, Beijing 100080, China 

e + e~ experiments producing charmonium are reviewed. It is found that the contribution of the continuum 
amplitude via virtual photon was neglected in almost all the experiments and the channels analyzed. It is shown 
that the contribution of the continuum part may affect the final results significantly in ip(2S) and -0(3770) decays, 
while the interference between continuum and resonance amplitudes may even affect the J/ip decays as well as 
the ip(2S) and ^(3770). This should be considered in analyzing the "p-zr puzzle" between J/ip and tp(2S) decays, 
and the difference between inclusive hadron and DD cross sections in 0;(377O) decays. 



1. Introduction 

There are three well-known problems in the 
study of the charmonium decays, namely the rel- 
ative phase between strong and electromagnetic 
amplitudes of the 1 charmonium decays, 11 pit 
puzzle" between J/ip and tp(2S) decays, and non- 
DD decays of 0(3770). 

The attempt to understand the strong decays 
of J/ip via three-gluon and the electromagnetic 
decays via one-photon annihilation reveals the 
relative phase between these two amplitudes is 
close to 90° I1I2I3I4| . while for the radially ex- 
cited ip(2S), the phase is 0° or 180° [ffi|. This 
indicates there would be no interference between 
these two amplitudes in J/ip decays, but strong 
interference in ip(2S) decays. 

It was found that in ip(2S) hadronic decays, 
some decay modes are abnormally suppressed 
compared with the corresponding J/ip decays 
based on perturbative QCD (pQCD) prediction. 
This suppression was first observed by the Mark- 
II in vector pseudoscalar (VP) decay modes like 
pir and K*K 5 , and confirmed by BES 0. 
Moreover, BES also observed the suppression in 
vector tensor (VT) decays of ip(2S) 0. This 
has led to active theoretical efforts in solving 
the problem |1I4I8I9| . Unfortunately, most of the 
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models were ruled out by the experiments, while 
some others need further experimental test. 

There is a renewed interest in 0(3770) stud- 
ies because of the upcoming high precision mea- 
surements by CLEO-c and BES-III (HJ. 
One of the puzzling problems in 0(3770) de- 
cays is that the DD cross section may be sig- 
nificantly lower than the inclusive hadronic cross 
section ^2]- This is in contradiction with the 
commonly accepted picture that -0(3770) decays 
predominantly to the OZI allowed DD states. 

These three topics play important roles in un- 
derstanding the charmonium decay dynamics. In 
this paper we examine what the experiments ob- 
serve and what theories analyze on charmonium 
produced in e + e~ experiments. We present a self- 
consistent analysis by considering the unavoid- 
able background process in e + e~ experiment, 
namely, the continuum process. We show that, 
for exclusive decays of these charmonium states, 
the contribution of this process could be very im- 
portant, or even if the direct contribution is rel- 
atively small, the interference between this term 
and other dominant amplitudes may contribute a 
non- negligible part. 

2. Experimentally observed cross section 

It is known that J/ip or ip(2S) decays into light 
hadrons via strong and electromagnetic interac- 
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tions. At the leading order in a s (m c ) and l/m c , 
it goes through three-gluon and one-photon an- 
nihilation of which the amplitudes are denoted 
by a-3g and a 7 respectively |2I13| . This is also 
true for -0(3770) in its OZI suppressed decay into 
light hadrons. In general, for the resonance IZ 
(K = J/tp, ip(2S) or -0(3770)), the cross section 
at the Born order is expressed as 



ctb(s) 



\a 3g 



(1) 



where yfs is the CM. energy, a is the fine struc- 
ture constant. If the J/0, ip(2S) or 0(3770) is 
produced in e + e~ collision, the process 
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hadrons 



(2) 



could produce the same final hadronic states as 
charmonium decays do ^3). We denote its am- 
plitude by a c , then the cross section becomes 



Ansa 



\a 3g 



(3) 



So what truly contribute to the experimentally 
measured cross section are three classes of di- 
agrams, i.e. the three-gluon decays, the one- 
photon decays, and the one-photon continuum 
process, as illustrated in Fig. ^ where the charm 
loops stand for the charmonium state, and the 
photons and gluons are highly off-shell and can 
be treated perturbatively. To analyze the exper- 
imental results, we must take into account three 
amplitudes and two relative phases. 

For an exclusive mode, a c can be expressed by 



a c (s) = ^e<* , 



(4) 



where 0' is the phase relative to a 3g ] ^F(s) de- 
pends on the individual mode, and for simplic- 
ity, the phase space factor is incorporated into 
\^F(s) | 2 . The one-photon annihilation amplitude 
can be written as 



a 7 (s) 



e v 



im n T t 



(5) 



where inn and T t are the mass and the total 
width of 7Z, T ee is the partial width to e + e~, is 




hadron 



e y (e e ) 



Figure 1. The three classes of diagrams of 
e + e~ — * light hadrons at charmonium reso- 
nance. The charmonium state is represented by 
a charm quark loop. 



the phase relative to 03^. The strong decay am- 
plitude a 3g is defined by C = \a 3 g/a 7 \, which is 



the relative strength to a 7 , so 



a 3g {s) = C 



ZT ee T(s)/(a.^/s) 
im-jiTt 



(6) 



For resonances, C can be taken as a constant. 

In principle, a 3gi a 7 and a c depend on individ- 
ual exclusive mode both in absolute values and 
in relative strengths. In this paper, for illustra- 
tive purpose, following assumptions are used for 
an exclusive hadronic mode: T(s) is replaced by 
y/R(s), where R(s) is the ratio of the inclusive 
hadronic cross section to the cross section 

measured at nearby energy |15| : in Eq. ffijl. 



C = 



' B(7Z — » ggg — ► hadrons) 
B(7Z — > 7* — > hadrons) 



(7) 



Here B(1Z 
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hadrons) = £? M + M -i?(s), 



where B^^- is the \i 



yT branching ratio; while 
B(7Z — > ggg — > hadrons) is calculated as fol- 
lowing: we first estimate the branching ratio of 
B(1Z — > jgg) + B(1Z — > ggg) by subtracting the 
lepton pairs, 7* — > hadrons, and the modes with 
charmonium production from the total branch- 
ing ratio (100%). Then using pQCD result pi] 
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Table 1 

Estimated amplitudes at J/tp, ip(2S) and -0(3770) peaks. 







m i>(2S) 


TO V(3770) 


l a 3 ff (m|)| 2 oc 
|a 7 (m^)| 2 oc 
\a c {m 2 n )\ 2 cx 


70% a J B fi ' 
13% a J J 4 ' 
20 nb 


19% 4 (25) 
1.6% 4 (2S) 
14 nb 


~ 1% 4 (3?70) 

2.5 x 10^ 5 ^ (3770) 
14 nb 



B(K -> jgg)/B(1Z -> .953) w 6% we obtain 
B(TZ — > ggg — > hadrons). Table ^ lists all the 
estimations used as inputs in the calculations, 
where er^ is the total resonance cross section of 



Born order at s = 



'TZ 



obtained from 



- m2 n) 2 + m n> 1 



1 p2 
TZ * 



(8) 



The cross section by e + e~ collision incorporat- 
ing radiative correction on the Born order is ex- 
pressed by |H| 



.00 



dxF{x, s)- 



o- (s(l - x)) 

1 -n(s(i 



(9) 



where er is ob or a' B by Eq. Q or 10, F(x,s) 
has been calculated in Ref. ^2] and n(s) is the 
vacuum polarization factor |18j : the upper limit 
of the integration x m = 1 — s m /s where ^/s^ is 
the experimentally required minimum invariant 
mass of the final state / after losing energy to 
multi-photon emission. In this paper, we assume 
that y/s^ equals to 90% of the resonance mass, 
Xffi — 0.2. 

For narrow resonances like J/tp and tp(2S), 
one should consider the energy spread function 
of e + e~ colliders: 
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2ttA 



2A2 



(10) 



where A describes the CM. energy spread of the 
accelerator, y/s and -s/s 7 are the nominal and ac- 
tual C. M. energy respectively. Then the experi- 
mentally measured cross section 



(11) 



CTexp{s) = / <7 r .c.{s)G{y/s, \fs l )dVs l 



49% and 29% respectively. The energy spread 
further reduces the cross sections of J/tp and 
tp(2S) by an order of magnitude. The radia- 
tive correction and energy spread also shift the 
maximum height of the resonance peak to above 
the resonance mass. Take tp(2S) as an exam- 
ple, from Eq. ©, a% (2S) = 7887 nb at tp{2S) 
mass; substitute ao(s) in Eq. © by ct^(s) in 
Eq. ijHJ), oy.c. reaches the maximum of 4046 nb 
at i/s = ?ti^(2s) + 9 keV; with the energy 
spread A = 1.3 MeV at BES/BEPC, combin- 
ing Eqs. (f51 — II 1|) . dexp reaches the maximum of 
640 nb at y/s = m^s) + 0-14 MeV. Similarly, 
at J/tp, with BES/BEPC energy spread A = 
1.0 MeV, the maximum of a exp is 2988 nb. At 
DORIS, the maximum of a exp at J/tp is 2190 nb 
(A= 1.4 MeV), and at tp{2S), it is 442 nb (A= 
2.0 MeV). In this paper, we calculate a exp at 
the energies which yield the maximum inclusive 
hadronic cross sections. 

To measure an exclusive mode in e + e~ exper- 
iment, the contribution of the continuum part 
should be subtracted from the experimentally 
measured cr' exp to get the physical quantity a exp , 
where <r exp and a' exp indicate the experimental 
cross sections calculated from Eas.lj§l-- tTT|) with 
the substitution of gb and a' B from Eqs. Q and 
respectively for uq in Eq. ©. Up to now, 
most of the measurements did not include this 
contribution and cr' exp = o exp is assumed at least 
at J/tp and 0(25). As a consequence, the theo- 
retical analyses are based on a exp , while the ex- 
periments actually measure a' 



exp " 



We display the effect from the continuum am- 
plitude and corresponding phase for J/tp, ip(2S) 
and V>(3770) respectively. To do this, we calculate 



The radiative correction reduces the maximum 
cross sections of J/tp, ip(2S) and 0(3770) by 52%, 
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Figure 2. (a) k as a function of <fi and </)' for 
ip(2S), with input from Tabled and (b) k as a 
function of <f> (<fr = <fi') for different ratios of |a3 9 | 
to |a 7 |: line 1 to 5 for a3 g = 0, |a3 g | = |a 7 |, 
|a 3s | = 3.4|a 7 |, |a 3s | = 5|a 7 | and \a 3g \ = 10|a 7 |, 
respectively. 



the ratio 



as a function of and </)', as shown in Fig. [2K 
for ip(2S) at Vs = ™^(2S) + 0.14 MeV for A = 
1.3 MeV. It can be seen that for certain values 
of the two phases, k deviates from 0, or equiva- 
lently the ratio o~' exp / 'cr exp deviates from 1, which 
demonstrates that the continuum amplitude is 
non-negligible. By assuming there is no extra 
phase between a 7 and a c (i.e. set <f> — 0'), we 
also work out the k values for different ratios of 
|o3 3 | to |et 7 |, as shown in Fig. \Bp: line 3 corre- 
sponds to the numbers listed in Table ^ line 1 
is for pure electromagnetic decay channels, and 
others are chosen to cover the other possibilities 



of the ratio |a 3ff | to |a~|. 

3. Continuum contribution for charmo- 
nium decay 

We now discuss separately the effect of contin- 
uum amplitude for 0(3770), ip(2S) and J/V>. 

At ■0(3770), the maximum resonance cross sec- 
tion of inclusive hadrons is 8 nb which predom- 
inantly decays into DD, while the continuum 
cross section is 14 nb which mainly goes to light 
hadrons. Assuming 1% of 0(3770) decays to non- 
DD interferes with the continuum amplitude, it 
could bring an effect of maximum 1.9 nb in the 
observed cross section. Such large constructive 
interferences could be responsible for the larger 
cross section of inclusive hadrons by direct mea- 
surement of e + e~ — > ■0(3770) — > hadrons than 
the DD cross section |12) . As to the exclusive 
decays, it could make some of the decay modes 
with small branching ratios more observable at 
the resonance. For example, if £5(0(3770) — > 
pir) ss 4 x 10 -4 (or equivalently, 0-^,(3770)^^ ~ 
0.003 nb) as suggested in Ref. and cr(e + e~ — > 
pir) w 0.014 nb at Born order by the model of 
Ref. |19j. then the maximum interference could 
be 0.011 nb, much larger than the pure contribu- 
tion from t/;(3770) decays. 

For 0(25), as can be seen in Fig. the ra- 
tio cr' exp /a exp could deviate from 1 substantially. 
In general, a^ g , a 7 and a c are different for differ- 
ent exclusive mode, so k could be different. This 
must be taken into account in the fitting of a 7 , 
a^g and the phase in between. It is noticeable 
that the observed cross sections of some electro- 
magnetic processes, such as tp(2S) — > 7r + 7r~, ujtt , 
and the famous puzzling ip(2S) — > pw, are three 
to four orders of magnitude smaller than the total 
hadronic cross section of the continuum process, 
which is about 14 nb. Form factor estimation |2T)] 
gives these cross sections at continuum compara- 
ble to the ones measured at the resonance |21| . 
It implies that a substantial part of the experi- 
mentally measured cross section could come from 
the continuum amplitude a c instead of the ip(2S) 
decays, and interference between these two am- 
plitudes may even affect the measured quantity 
further. Therefore it is essential to measure the 
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production rate of 7r + 7r~, W7r° and pir at the con- 
tinuum in order to get the correct branching ra- 
tios of the ip(2S) decays. The same holds for VT 
decays of ip(2S). 

As for J/tp, the interference between the am- 
plitude a c and the resonance is at the order of 
a few percent. It is smaller than the statistical 
and systematic uncertainties of current measure- 
ments. Nevertheless, for future high precision ex- 
periments such as CLEO-c TO] and BES-III [TT] . 
when the accuracy reaches a few per mille or even 
smaller level, it should be taken into account. 

4. Dependence on experimental conditions 

Here we emphasize the dependence of the ob- 
served cross section in e + e~ collision on the ex- 
perimental conditions. The most crucial ones are 
the accelerator energy spread and the beam en- 
ergy setting for the narrow resonances like J/ip 
and ip(2S). 

Fig. |2| depicts the observed cross sections of 
inclusive hadrons and pairs at ip(2S) in 

actual experiments. Two arrows in the figure 
denote the different positions of the maximum 
heights of the cross sections. The height is re- 
duced and the position of the peak is shifted due 
to the radiative correction and the energy spread 
of the collider. However, the energy smear hardly 
affects the continuum part of the cross section. 
The p + fi~ channel is further affected by the in- 
terference between resonance and continuum am- 
plitudes. As a consequence, the relative contribu- 
tion of the resonance and the continuum varies as 
the energy changes. In actual experiments, data 
are naturally taken at the energy which yields the 
maximum inclusive hadronic cross section. This 
energy does not coincide with the maximum cross 
section of each exclusive mode. So it is important 
to know the beam spread and beam energy pre- 
cisely, which are needed in the delicate task to 
subtract the contribution from a c . 

It is worth noting that in principle if a c is not 
considered correctly, different experiments will 
give different results for the same quantity, like 
the exclusive branching ratio of the resonance, 
due to the dependence on beam energy spread 
and beam energy setting. The results will also be 
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Figure 3. Cross sections in the vicinity of ip(2S) 
for inclusive hadrons (a) and /i + pT (b) final 
states. The solid line with arrow indicates the 
peak position and the dashed line with arrow the 
position of the other peak. In (b) , dashed line for 
QED continuum (<r c ), dotted line for resonance 
(<r R ), dash dotted line for interference (a 1 ), and 
solid line for total cross section(er Tot ). 



different for different kinds of experiments, such 
as production of J ftp and ip{2S) in pp annihila- 
tion, or in B meson decays. This is especially 
important since the beam spreads of different ac- 
celerators are much different [2] and charmonium 
results are expected from B-factories. 

5. Summary and perspective 

In summary, the continuum amplitude a c , by 
itself or through interference with the resonance, 
could contribute significantly to the observed 
cross sections in e + e~ experiments on charmo- 
nium physics. Its treatment depends sensitively 
on the experimental details, which has not been 
fully addressed in both e + e~ experiments and 
theoretical analyses. So far, most of the mea- 
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surements have large statistical and systematic 
uncertainties, so this problem has been outside 
the purview of concern. Now with large J/ip and 
ip(2S) samples from BES-II [22] and forthcom- 
ing high precision experiments CLEO-c [TO] and 
BES-III [TT], the effect of a c needs to be treated 
properly. To study it, the most promising way 
is to do energy scan for every exclusive mode in 
the vicinity of the resonance, so that both the 
amplitudes and the relative phases could be fit 
simultaneously. In case this is not practicable, 
data sample off the resonance with comparable 
integrated luminosity as on the resonance should 
be collected to measure |a c |, which could give an 
estimation of its contribution to the decay modes 
studied. The theoretical analyses based on cur- 
rent available e + e" data, particularly on il>(2S) 
may need to be revised correspondingly. 
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